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P.5 A NUMERICAL STUDY OF SCALAR GRADIENTS IN 
KELVIN-HELMHOLTZ BILLOWS 

J. W. Parker 

Department of Electrical and Computer Engineering 
University of Illinois, 1406 W. Green, Urbana, IL 61801 

S. A. Bowhill 

Department of Electrical Engineering 
University of Lowell, 1 University Avenue, Lowell, MA 01854 

A high-resolution numerical technique is used to model the development of a 
periodically perturbed shear layer imbedded in an initially vertical gradient of a passive scalar. 
The technique follows the development of the vorticity through an initial linear growth stage 
and well into the nonlinear development of Kelvin-Helmholtz billows, in the zero-viscosity, 
zero-diffusion limit. The resulting scalar distribution rapidly develops regions of extremely 
sharp scalar gradients, which wind around the periodically spaced vortical low-gradient cores. 
Vertical cross sections through different parts of the billow structure are presented and 
compared with rocket measurements of electron density fine structure in the mesosphere. 
Gradient limits imposed by finite diffusion are calculated, and implications for atmospheric 
radar observations are discussed. 


What is the Kelym-Mel mholto Instability 2 


A region of fluid shear is unstable. 


Disturbances grow, billows form. 



Rockets, radar observe advected scalar quantities (temperature, ionization). 
What happens to scalar distribution? 


What features should show up in rocket, radar data? 
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Figure 1. Billow amplitudes vs. time scaled by growth factor for kh = 0.15, 0.25, 0.35, 0.45, 
055, 0.65, 0.75, 0.85, 0.95, 1.05, 1.15, and 1.25. 
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Figure 2. Scalar development in billow with time, (a) initial condition, with constant scalar 
vertical gradient, contours at yfk = ± 0.06 are also boundaries of vorticity region, (b) tea = 
14.0. (c) tco = 19.6. (d) tto = 25.2. (e) tco = 30.8. (f)t© = 36.4. 
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Figure 3. Vertical cross sections of scalar for pattern of Figure 2(f), through stagnation point, 
two intermediate locations between stagnation point and core, and core center. 
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Figure 4. Log average of six energy spectral densities from cuts through scalar irregularities in 
simulated billows. 
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Figure 5. Reconstruction of digital electron density data for a portion of 28 May 1975 Peru 
rocket shot [Stoltzfus et al., Adv. Space Res., 4(6), 143, 1984]. 



Figure 6. Power spectrum from 28 May 1975 Peru rocket shot. Spectrum is average of 
spectra from two parts of data of Figure 4, with noise estimate subtracted, and smoothing in 
wave number domain. 





